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Abstract 

We consider the ultradiscretization of a solvable one-dimensional chaotic map which arises from the duplication 
formula of the elliptic functions. It is shown that ultradiscrete limit of the map and its solution yield the tent map 
and its solution simultaneously. A geometric interpretation of the dynamics of the tent map is given in terms of the 
tropical Jacobian of a certain tropical curve. Generalization to the maps corresponding to the ra-th multiplication 
formula of the elliptic functions is also discussed. 



1 Introduction 

In this article, we consider the following map 

f . , 4Z„(1 - Zn)(l ~ khn) n n 

Zn+1 = f(Z„) = (l_fc2 z 2)2 ' (L1) 

which admits the general solution 

z n = sn 2 (2" Mo ;fc), (1.2) 

describing the orbit in [0, 1]. Here m(u;k) is Jacobi's sn function, < k < 1 is the modulus, and uq is an arbitrary 
constant. In fact, (11. U can be reduced to the duplication formula of sn function: 

2sn(«; k) cn(w; k) dn(«; k) 

sn(2u;k) = , (1.3) 

1 - k z sn*(u; k) 

cn 2 (u;k) = l-sn 2 (u;k), dn 2 (u;k) = 1 - k 2 sn 2 (u;k), (1.4) 

where cn(«; k) and dn(t<; k) are Jacobi's cn and dn functions, respectively. The map ( 11.11 ) is a generalization of the 
logistic map (or Ulam-von Neumann map): 

Zn+i = 4z„(l - z n ), Zn = sin 2 (2 n M ). (1.5) 

The map (11. Il l was first considered by Schroder ||28l in 1871, and it has been studied by many authors [10, 12 38l [39l . 
It is now classified as one of the (flexible) Lattes maps ll2D . In this article, we call ( 11. It the Schroder map. 
It is well-known that the Schroder map is conjugate to the tent map for X„ e [0, 1] 



X„+\ = T2(X n ) =1-2 

Namely, we have the relation 



1 

Xn~- 



\2X„ 

12(1 -X n ) ±<X n <l. 



sofo S - l =T 2 , s(z) = -J—sn-\y/i;k), (1.7) 
K(k) 

where K(k) is the complete elliptic integral of the first kind 

f 1 dx 

K(k)= — . (1.8) 

Jo V(! -x 2 )(1 -k 2 x 2 ) 

The purpose of this article is to establish a new relationship between the Schroder map and the tent map through 
a certain limiting procedure called the ultradiscretization J35|. The method of ultradiscretization has achieved a great 
success in the theory of integrable systems. From the integrable difference equations, various interesting piecewise 



1 



linear dynamical systems have been constructed systematically, such as the soli ton cellular automata lPfll6l [T7ll22"l|29l 
[30] [31] [34J [371 ST] and piecewise linear version of the Quispel-Roberts-Thompson (QRT) maps l|23l |26l [32] [36l . The 
resulting piecewise linear discrete dynamical systems can be expressible in terms of the max and + operations, which 
we call the ultradiscrete systems. The key of the method is that one can obtain not only the equations but also their 
solutions simultaneously. It also allows us to understand the underlying mathematical structures of the ultradiscrete 
systems g] E \U\ [TJ M CES M 121] S3 . 

In this article, we apply the ultradiscretization to the Schroder map dl.lt and its elliptic solution jl.2\ . As a result 
they are reduced to the tent map and its solution. We also clarify the tropical geometric nature of the tent map; we 
show that the tent map can be regarded as the duplication map on the Jacobian of a certain tropical curve. 

2 Ultradiscretization of the Schroder map 

The key of the ultradiscretization is the following formula: 

/ A £ \ 

lim elog <?' + e e H = max(A, B, ■ ■ •), (2.1) 

where the terms in log must be positive, and the dominant term survives under the limit. We note that the orbit of the 
map (11. U is always restricted in [0, 1] if the initial value is in this interval. Since this is somewhat too restrictive for 
ultradiscretization, we apply the fractional linear transformation 



Zn 



Zn 



1 -Zb' 



which maps [0, 1] — > [0, oo). Then the Schroder map Jl.lt and its solution ( 11.2b are rewritten as 

4jc„(1 + x n )(l +£' 2 x„) 



X n +l = 4>( x n) = 



sn 2 (2" Mo ; k) sn 2 (2" Mo ; k) 



l-Zn 1 - sn 2 (2"M ; k) cn 2 (2"« ; k) 

respectively. We note that the map (12.31 l can be obtained from ( II. Il l by replacing as z n 
On the level of solution, this corresponds to Jacobi's imaginary transformation 



— i sn(iu; k') = 



sn(u; k) 
cn(u; k) 



(2.2) 

(2.3) 
(2.4) 

k' = VT^P. 
(2.5) 



Figure Q] shows the map functions of ( II. lb and ( 12.31 ). Note that f(z) and <p(x) have poles at z = +l/k and x = +l/k', 
respectively. 




1 2 3 4 5 



Figure 1: Map functions of CD) (left: k = 0.7) and (Bright: k' = 0.8) 



Now we put 



x„ = exp 



, k' = exp 



L 
Ye 



, (0<Jf < 1, L>0). 



(2.6) 



Then ( 12. 3t is rewritten as 



X n+l = F e (X n ) = elog 



(2.7) 



Taking the limit e — > +0 by using the formula ( 12. U . we obtain 

X„ +1 = F(X n ) =X„+ max(0, X„) + max(0, X„ - L) - 2 max(0, 2X„ - L) 

'X n X n < 0, 

2X„ < X„ < §, 

-2X„ + 2L | < X„ < L, 

-X„ + L L < Z„. 



(2.8) 




Figure 2: Left: map function of the ultradiscrete Schroder map (12.81 1. Right: limit transition of the map function F e {X) 
for L — 1.5. Dashed line: e = 0.3, dot-dashed line: e = 0.1, solid line: e = 0.01. 



Remark 2.1 Although the terms in log m the formula ( 12.71 ) mi«f fee positive in general, the negative terms can also 
exist as long as they are not dominant in the limit. For example, we have 

I A £\2 / 2/1 A+B 2B\ 

lim elogle* -e e ) = run elog [e * -2c « + e « = 2 max(A, B). (2.9) 

v ' e-H-0 ^ ' 

We call the map (12 . 8b the ultradiscrete Schroder map. Figure [2] shows the map function of d2.8l l and limit transition 
of the function F e (X). The dynamics of the map d2.8l l is described as follows: if the initial value Xo is in [0, L], the 
map is the tent map and X„ e [0, L] for all n. If Xo e (— oo, 0], then X„ = Xo for all n > 1. Finally if Xq g [L, oo), then 
X\ = -Xo + L < and X„ = X\ for all n > 1 . Therefore the ultradiscrete Schroder map ( 12.81 ) is essentially the tent map 
on [0, L] 

X„_l 

L 2 



X„+i - I I - 2 



X n €[0,L], (2.10) 



and otherwise the dynamics is trivial. 

Now let us consider the limit of the solution by using the ultradiscretization of the elliptic theta functions 113211 (see 
also Ifl4ll24ll25ll ). Jacobi's elliptic functions are expressed in terms of the elliptic theta functions #,(v) (i = 0, 1, 2, 3) 

sn(u;k)= , cn(u;k) , (2.11) 

#2(0)# (v) # 2 (0)W 

u=«m0))\ * = (^)\ (2-12) 
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where 



if (v) - 
0i (v) = 


\ 1 / i \n n 2 Jin 
" 2-A ^ q Z ' 

neZ 

neZ 


(2.13) 
(2.14) 


0a(v) = 
3 (v) = 


neZ 


(2.15) 
(2.16) 



and z = exp[inv]. We parametrize the nome q as 

q = exp 



67T 



, 0>O. 



(2.17) 



Applying Jacobi's imaginary transformation (or Poisson's summation formula) the elliptic theta functions are rewritten 

as 

„ ■ ■ ■ ' 2 

^2> P 



My) = 

01 (V) = 

2 (v) = 

#3(v) = 



neZ 

-£(-!)" exp 



neZ 



— (y-»X 
e 



X exp - 

neZ 



■ (v-«y 



Asymptotic behaviour of these functions for e — > +0 is given by 



00(0) 

02(0) 
03(0) 

(0o(v)) 2 
(01 (v)) 2 



e r 

— exp - — 

67T L 4e 

VI( i+2ex pH])' 



— exp 

en 

8 

— exp 
en 

6 





((v)) - 






ill 












-II 


((v)) - 









— {((v))r 

L e 



- exp 







— {((v»-ir 

L e 



(0 2 (v)) 2 ~ —(exp 
en \ 

where ((v)) is the decimal part of v, namely, 

((v)) = v - Floor(v), < ((v)) < 1. 

Then we have 

L~\ 1 ,, 2 t /0 2 (O)\ 4 16exp[-f](l + 4exp[-f1) 



k' = exp 



x„ = exp 



U(0)j (l + 2exp[-f]) 4 
X„ sn^ M ;fc) /0 3 (O)0!(v)\ 2 (l + 2exp [-f ]f exp [»] 



cn2(«;fe) \0o(O)0 2 (v)/ 4 (l - exp [f [2((v)) - l]]) 2 
4 



(2.18) 
(2.19) 
(2.20) 
(2.21) 

(2.22) 
(2.23) 
(2.24) 
(2.25) 

(2.26) 
(2.27) 

(2.28) 



which yield in the limit e — > +0 



L = 

X„ = 







h 


«v))-i) 



2"vo, 



(2.29) 
(2.30) 



respectively, where vo is an arbitrary constant. We note that in taking the limit of x n , we have put the arbitrary constant 

wo as 

6 

m = -v (2.31) 
e 

so that 

2"u 



= 2"v 



2"v (e-»+0). 



(2.32) 



^■(i? 3 (0)) 2 7!-(#3(0)) 2 

One can verify that ( 12.291 ) and ( 12.301 1 actually satisfy the ultradiscrete Schroder map ( 12.8b or ( 12.101 ) by direct calcula- 
tion. Therefore we have shown that through the ultradiscretization the Schroder map ( 12.31 ) and its solution ( 12.4b yield 
the map ( 12.8b (or ( 12.10b ) and its solution (12.301 ) simultaneously. 



Remark 2.2 



(1) The fundamental periods of ^"j!^ & re 2K(k) and 2iK(k'). In the ultradiscretization of the elliptic theta functions, 
we have parametrized the name q as (12.17\) , which implies that the ratio of half-period T is given by r = i~ ant/ 



*(*) = ^(# 3 (0)) 2 ~ K(k') = -^(# 3 (0)) 2 r = ^(# 3 (0)) 2 ~ |, 



2e 



26» 



(2.33) 



ase-> +0. S/nce we have u = -v, the fundamental periods with respect to v tend to 1 and j(isf-» +0. This 
implies that the ultradiscretization of the elliptic functions is realized by collapsing the imaginary period and 
keeping the real period finite. 



(2) The Schroder map ( 17.71 ) is reduced to the logistic map ( 17.51 ) for k — 0. This corresponds to the ultradiscrete 
Schroder map ( 12.81 ) with L = 0, 

X„ +l = -\X n \, (2.34) 

whose dynamics is trivial, and the solution ( 12.301 ) becomes X n — 0. Therefore ultradiscretization of the logistic 
map does not yield an interesting map [9 J. In fact, we see that this case is not consistent with the ultradiscrete 
limit, since the asymptotic behaviour ofK(k) and K(k') as k — > is given by 



K(k) ~ ~ K(k') ~ log 1 

2 K 



(2.35) 



One can apply the same procedure to the following map which originates from the triplication formula of sn 2 lfl2l 
EDE11 

z n [k A z A n - 6k 2 z 2 n + 4(k 2 + \)z„ - 3} 2 



Zn+l = g(Zn) 



which is rewritten as 



x„+i = y(x n ) = 



{3k 4 zi - 4k 2 (k 2 + l)z 3 „+ 6k 2 z 2 - 1} 

x n {k' 4 x 4 n - 6k' 2 x 2 n - 4(k' 2 + Y)x n - 3} 
[ik' A 4 + 4k' 2 (k' 2 + \)xl + 6k' 2 x 2 n - 1 



z„ = sn (3"m ; k), 



sn-(3"M ; k) 
cn 2 (3"M ;/t)' 



(2.36) 



(2.37) 



by the transformation (12.21 ). The map functions g(z) and y(x) are illustrated in figure [3] Then ultradiscretization of 
(12.37b yields the map 

X n+1 = G(X„) = X„+2 max(0, X n , 4X„ - 2L) - 2 max(0, 3X„ - L, 4X„ - 2L) 
'x„ X„ < 0, 

3X„ 0<X„< |, 

-3X„ + 2L ^<X n <^, 



3X n - 2L 
X„ 



3 
2L 



(2.38) 



< X„ < L, 



L < X n , 
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Figure 3: Map functions of ( 12361 1 (left: k = 0.7) and fl237]>(right: k! = 0.8) 



and its solution 









«v»-i) 



I, v = 3"v . 



Figure |4] shows the map function G(X„) and the limit transition of the map function of 



X n +\ = G e (X n ) = elog 



-Vh I Mn-IL , 2X tl -L 



-6e< -4(e-- + l)ef -3 



(2.39) 



(2.40) 



{3e^~ +4(e~~ +e~-))e— + 6e^" - 1} 
We note that one can directly ultradiscretize the map ( 12.36b to obtain ( 12.381 ). however, the solution x„ = sn 2 (3 n Mo; k) 





Figure 4: Left: map function of the map ( 12.38b . Right: limit transition of the map function G e (X) for L — 1.5. Dashed 
line: e - 0.3, dot-dashed line: e — 0.1, solid line: e = 0.01. 



degenerates to the trivial solution X„ = 0. Thus it is important to consider ( 12.371 ) in order to obtain the limit which is 
consistent with the solution. 

It is possible to apply ultradiscretization to the maps arising from the m-th multiplication formula of sn 2 lfT2ll2Tll 
in a similar manner. 



3 Geometric description in terms of the tropical geometry 

It is shown in [5, 25 1 that the tropical geometry provides a geometric framework for the description of the ultradiscrete 
integrable systems. Therefore it may be natural to expect that a similar framework also works well for our case. In 
this section, we show that the ultradiscrete Schroder map can be interpreted as the duplication map on the Jacobian of 
a certain tropical curve. As for the basic notions of the tropical geometry, we refer to lfTll7l l27l . 
We first consider the elliptic curve 

[xy - b(x + y) + c] 2 = 4d 2 xy, (3.1) 
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parametrized by 



(x,y) = 



sn 2 (i<; k) sn 2 (« + rj\ k) 
cii 2 (m; k) ' cn 2 (t< + J];k)j' 



where rj is a constant and a, b, d are given by 

1 cn 2 (j];k) 1 d _ 1 dnfafc) 

~ k 12 sn 2 (?;; ifc) ' C_ fe ;2 ' ~ ~k^sn 2 (r];ky 

respectively. Eliminating rj in ( 13.31 ). we see that b and a! satisfy the relation 

k' 2 d 2 = (I + k' 2 b)(l + b). 



(3.2) 



(3.3) 



(3.4) 



We may regard the Schroder map ( 12.31 ) as the projection of the dynamics of the point on the elliptic curve ( 13. U to the 
x-axis. 

We next apply the ultradiscretization to the elliptic curve. Putting 



1 



— - — 2 — 

and taking the limit e -> +0, ( 13.11 ) and ( 13.41 ) yield 

max(2X + 2F, B + 2X, B + 2Y, 2L) = X + Y + D 

and 

-L + D = max ^0, — - Lj + max ^0, — j , 
respectively. The condition ( 13.7b gives the following three cases: 

(i) B > 2L > 0, D = B, 

(ii) 2L>B>0, D = L+-, 

2 



(iii) 



> B, D - L. 



(3.5) 

(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 



For each case, the set of points defined by (13.6b is (i) a line connecting (f » f ) and (L - j,L - f ), (ii) a rectangle 
with vertices (0, L - f ), (L - f , 0), (L, f ) and (f , L), (iii) a line connecting (f ,L - f ) and (L - f , |), respectively, as 
illustrated in figure|5] In the following, we consider only the case (ii) and we denote the rectangle as C. 






y 






B 




L 

N. L 


L-°- 
2 












B 

2~ 







Figure 5: Ultradiscretization of the elliptic curve ( 13.11 ). Left: case (i), center: case (ii), right: case (iii). 
Let us recall some notions of the tropical geometry. The tropical curve defined by the tropical polynomial 

E(X, Y) = max {A( au a z ) + a x X + a 2 Y), 31 e Z 2 , 



(3.11) 



is a set of points (X, Y) e R 2 where S is not smooth. Here 31 is a finite subset of Z 2 called the support, and we denote 
as A(3l) the convex hull of J{. Let be the triangle in Z 2 with vertices (0, 0), (d, 0), (0, d). Then the degree of the 
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tropical curve is d if A(Jl) is inside F</ but not inside T^-i [40|. The genus of the tropical curve is defined as the first 
Betti number of the curve, namely the number of its cycles IfTl [T8l [T9l . 
We consider the tropical polynomial 

x ¥(X,Y) = max(2X + 2Y,B + 2X,B + 2Y,2L,X+Y + D), (3.12) 

under the condition d3.9t . Let C be the tropical curve defined by T, which is illustrated in figure [6] Then the degree 
and the genus of C are 4 and 1, respectively. Note that the rectangle C is exactly the cycle of C. 




Figure 6: Left: tropical curve C defined by ( 13.12b . Right: support of ( 13.12b . 



Vigeland [40] has successfully introduced the group law on the tropical elliptic curve. Unfortunately, however, 
his definition of tropical elliptic curve is limited to "smooth" curve of degree 3 and hence it does not cover our case. 
Nevertheless, it is possible to define the tropical Jacobian J(C) of C 15] 1201 and characterize the dynamics of the 
ultradiscrete Schroder map ( 12.101 ) on it in the following manner: let V, and E, {i = 1 , . . . , 4) be the vertices and edges 
of C defined by 

Vi =O = (0,L-1),V 2 = (l-^,0),V 3 = (l,I),V 4 = (1,l), (3.13) 

ViV 2 = E u V 2 V 3 = E 2 , V 3 V 4 = E 3 , V 4 Vi=£ 4 , ( 3 - 14 ) 
respectively. The length of each edge is given as 

m = V2"(z, - |) , \E 2 \ = ^p, \E 3 \ = V2"(l - I) , \E A \ = ^B. (3.15) 

The primitive tangent vector for each edge is 

V! = (1,-1), v 2 = (l,T), v 3 =(-l,l), v 4 = (-l,-l). (3.16) 

We introduce the total lattice length X. as the sum of the length of each edge scaled by the length of corresponding 
primitive tangent vector, which is computed as 

£=T^=2L. (3.17) 
it l v ''l 

Then the tropical Jacobian J(C) is defined by 

J(C) = R/XZ = R/2LZ. (3.18) 
The Abel-Jacobi map [i : C — > J(C) is defined as the piecewise linear map which is linear on each edge satisfying 

Ai(Vi) = 0, n{V 2 ) = L~^> M^ 3 ) = L, fx(V 4 ) = 2L - |. (3.19) 

Let7r:C — > Rbe the projection of the point on C to the X-axis. Letpbe the map defined by p = nofT 1 ; J(C) — > R 
which maps fi(P) (P e C) to the X-coordinate of P. Here we note that n~ l is 1:2 and we define n~ x {X) to be the point 
on C whose F-coordinate is smaller. In this setting, p{p) (p e 7(C)) can be written as 

i [p < p < L, 

p(p) = (no^)(p)=\ I ~ '- ' (3.20) 

\-p + 2L L < p < 2L, 
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A 




Figure 7: Left: correspondence between X and 7(C) by p. Right: duplication map ip 2 and <&2- 



as shown in the left of figure [7] 

Now we define the duplication map ip 2 : 7(C) — > 7(C) by 

<p 2 (p) = 2p (mod £), p e 7(C), (3.21) 

and introduce <$>2 '■ R — > R as the conjugation map of tp 2 by p, 

3*2 = P ° V2 °P _1 - (3.22) 

In order to write down the map O2 explicitly, we introduce p',p" e 7(C) for P — (X, Y) e C by 

p' =p-\X) = (jxo 7 r- 1 )(X)=X, p" =<p 2 (p') = 2p' =2X. (3.23) 

Then the map $2 is expressed as follows (the right of figure|7]): 

(1) For < X < \ : since < p" < L, ( f3T20b implies 

<t> 2 (X)=p(p") = 2X. (3.24) 

(2) For § < X < L: since L < p" < 2L, ( T3~20l i implies 

G>2(X) = p(p") = -2X + 2L. (3.25) 

The dynamical system 

Xn_ 1 

L 2 



Z n+1 = <1>2(X„) = L 1 - 2 



2X„ 0<X<§, 
-2X„ + 2L k <X < L, 



(3.26) 



coinsides with the ultradiscrete Schroder map ((2T0]l. Therefore we have shown that the ultradiscrete Schroder map 
( 12.1 Oi l can be regarded as the duplication map on the Jacobian 7(C) of the tropical curve C defined by the tropical 
polynomial ( 13.12b . 

Similarly, we define the triplication map tpj : 7(C) — > 7(C) by 

Wip) = 3 P (mod £), p e 7(C), (3.27) 

and introduce (D3 : R — > R as the conjugation map of ^3 by p, 

<E> 3 =po(p 3 op- 1 . (3.28) 

Then the corresponding dynamical system is given by 

\3X n 0<X„< £, 



X n+1 = <S 3 (X„) = 



-3X„ + 2L §<X<f, 
3X„ - 2L ^ <X<L 



3 

= 3X„ - 2 max(0, 3X„ - L) + 2 max(0, 3X„ - 2L), (3.29) 
which is equivalent to ( 12.38b on [0, L]. For general to, the TO-th multiplication map yields the dynamical system 

m-l 

X„ + , = <b m (X n ) = toX„ + 2 ^(-1)'* max(0, mX„ - iL), (3.30) 

;=1 

which may be regarded as the ultradiscretization of the map arising from the TO-th multiplication formula of ^ . 
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4 Concluding Remarks 

In this article, we have presented a new relationship between two typical chaotic one-dimensional maps, the Schroder 
map and the tent map, through the ultradiscretization. Although the ultradiscretization has been developed in the 
theory of integrable systems, the results in this article imply the possibility of applying the method to wider class of 
dynamical systems. Our results also suggest that the tropical geometry combined with the ultradiscretization provides 
a powerful tool to study a piecewise linear map, since the ultradiscretization translates the geometric background of the 
original rational map into that of the corresponding piecewise linear map. It would be an interesting problem to study 
various ultradiscrete or piecewise linear systems, such as ultradiscrete analogues of Painleve systems, generalized 
QRT maps, and higher-dimensional solvable chaotic maps in this direction. 
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